Enhancement of GMR due to spin-mixing in magnetic multilayers with a 

superconducting contact 
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We study the Giant Magnetoresistance (GMR) ratio in magnetic multiayers with a single super- 
conducting contact in the presence of spin-mixing processes. It has been recently shown |jj that 
the GMR ratio of magnetic multilayers is strongly suppressed by the presence of a superconducting 
contact when spin-flipping is not allowed. In this Letter we demonstrate that the GMR ratio can 
be dramatically enhanced by spin-orbit interaction and/or non-collinear magnetic moments. The 
system is described using a tight-binding model with either s-p-d or s-d atomic orbitals per site. 



PACS numbers: 75.70.Pa, 74.80.Dm 

Hybrid nanostructures form a fascinating melting pot 
for studying the interplay between fundamental quantum 
phenomena, often revealing new and unexpected physics. 
One recently-recognized class of such structures, involv- 
ing the coexistence of superconducting contacts and fer- 
romagnetic domains, has led to the identification of a 
number of fundamental issues [|]-^) > several of which are 
currently unresolved. In this Letter, we examine one 
such issue, posed by experiments on giant magnetoresis- 
tance (GMR) in magnetic (M) multilayers with supercon- 
ducting (S) contacts and current-perpendicular-to-the- 
plane (CPP). Recognizing that the sub-gap conductance 
of such structures is mediated by Andreev scattering, it 
was recently noted jl| that in the absence of spin-flip 
processes, the conductance of a metallic (i.e. diffusive) 
multilayer in the presence of aligned magnetic moments 
is almost identical to that of the multilayer when adja- 
cent moments are anti-aligned and therefore the conven- 
tional GMR ratio should be strongly suppressed. Since 
large GMR ratios are observed experimentally ||, it is 
clear that even a qualitative understanding of transport 
in such structures must incorporate the effects of spin- 
mixing. The aim of this Letter is to present the first 
theoretical description of CPP GMR in M-multilayers 
with S-contacts, which incorporate spin-flip scattering. 
As sources of spin-mixing we consider both spin-orbit 
(SO) coupling and non-collinear magnetizations in adja- 
cent magnetic layers. 

The system under consideration is a disordered mag- 
netic multilayer consisting of an alternating sequence of 



magnetic layers each of length Im and non-magnetic lay- 
ers (N) of length Im- The building block of the magnetic 
structure is the bilayer [M/N] of length Is — In + Im- 
The magnetic moments of even-numbered M-layers make 
an angle 9 relative to those of odd-numbered M-layers. 
Experimentally, 9 can be varied by applying an external 
magnetic field with antiparallel (AP) alignment (9 = tt) 
typically occurring at zero field and parallel (P) align- 
ment (9 = 0) at large enough fields. The current flows 
perpendicular to the planes of the multilayer, which 
makes contact with a metallic normal lead on the left- 
hand side of the multilayer and a superconducting lead 
on the right-hand side. GMR is the drastic increase in 
electrical conductance G(9) that occurs when the system 
switches from the AP to the P alignment with the con- 
ventional GMR ratio defined by: p = ^^q^§^-- 

Following jl) , the multilayer and leads are modelled us- 
ing a tight-binding Hamiltonian on a cubic lattice with 
hoppings to nearest neighbours. Lattice imperfections 
and impurities are simulated by adding to the on-site en- 
ergies a random number in the range [— ttj+tt]. The 
on-site Hamiltonian has the following structure: 
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where ffPT(i) [ s the Hamiltonian for up (down)-spin par- 
ticles (s and d bands), H h ^ = -iPKD* [ s the Hamil- 
tonian for up(down)-spin holes and A is the supercon- 
ducting order parameter. Here [i xy = —jjL x + ijJ, v , where 
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Hx(y) is the a;(y)-component of the exchange field fx. Note 
that pi is non-zero only for electrons in the d-band in the 
M-layers and A is non-zero only in the right-hand-side 
superconducting lead. Within the tight-binding formu- 
lation, SO interaction can be included by adding to the 
Hamiltonian the following term: 



V so = Vi 



R ■r crt r~ CT 



(2) 



where V\ is a constant which determines the interaction 
strength, a is a vector of Pauli matrices, Rij is the unit 
vector which connects site i with the neighbouring site j. 
c a a i is the annihilation operator for electrons of spin a in 
the a (s, d)-band on site i. In the presence of disorder 
(|J) produces spin-flip scattering since it couples electrons 
with different spin on neighbouring sites. 

In the presence of disorder, to study the largest pos- 
sible sample cross sections, we consider 2 orbitals per 
site, which is the minimal model capable of reproducing 
scattering potential at the N/M interface and interband 
scattering |[o) . The tight-binding parameters are chosen 
to reproduce the GMR ratio and conductances obtained 
from an ab initio material specific calculation for Cu/Co 
multilayers Jll| . 

In the presence of spin-flip scattering the two-spin-fluid 
approximation does not hold and the zero-temperature, 
zero-bias, normal-state Landauer formula takes the form: 



qNN _ 



^E Tr { r ' tr '} 



(3) 



where t aa is the matrix of transmission amplitudes for 
injected ex'-spin electrons in the left-hand lead into cr-spin 
electrons in the right-hand lead. When the right-hand 
lead is in the superconducting state, the conductance is 
given by p) 



G 



NS 



(4) 



where r° a is the Andreev reflection matrix for injected 
cr'-spin electrons in the left-hand lead to be reflected 
into cr-spin holes. In what follows, the scattering ampli- 
tudes are calculated exactly by solving the Bogoliubov- 
de Gennes equation using an efficient recursive Green's 
function technique EUllJjl^]. 

We shall now turn to the central results of this Let- 
ter, namely that in the presence of strong-enough spin- 
mixing, either produced by SO coupling or non-collincar 
moments, GMR in the presence of a S-contact approaches 
that of two normal contacts, with values of p of the order 
of 100 %. First consider the effect of SO coupling within 
the multilayer. Fig. [l] shows the conventional GMR ratio 
as a function of the SO interaction strength Vy for a disor- 
dered multilayer of 40 bilayers, with Im = 15 and In = 8. 



As expected, in the NN case p decreases monotonically 
from ~ 200 % at Vi = 0, to zero at large Vi (~ 0.17 eV). 
In contrast for the NS case, p initially increases with in- 
creasing Vi, eventually joining the NN curve at Vi ~ 0.08 
cV. As a second source of spin-mixing, consider the ef- 
fect of non-collinear magnetic moments when Vy = 0. 
Fig. shows the ^-dependence of the GMR ratio defined 
as p{6) = G{ %^ ] ■ Whereas in the NN case p(9) de- 
creases monotonically with increasing 6, in the NS case 
p(9) exhibits a pronounced maximum around 9 = n/8. 



To understand these results, first consider the case of 
non-collinear moments. In the presence of two normal- 
metallic contacts, the conductance G(9) has been the- 
oretically studied in [fli 15 where it is predicted that 



G(9) — G(ir) tends monotonically to zero as 9 varies from 
to 7T. In addition, the dependence of the resistance 
on the angle 9 has been experimentally found |16| to 
contain a term proportional to cos 2 (9/2) and a second 
term proportional to cos 4 (9/2). In the presence of a S- 
contact, where G(0) ~ G(ir), this behaviour is drastically 
changed by the presence of an extremum which occurs 
at some intermediate angle 9 C , the value of which de- 
pends on the interplay between competing effects. Since 
9(H) is a function of the applied magnetic field H, the 
presence of a S-contact introduces a new characteristic 
field H c for which 9(H C ) = 9 C . For a disordered multi- 
layer of 22 bilayers, with Im — 30, In = 16, the insert 
in Fig. U shows the ^-dependence of the conductance 
divided by the number of open channels in the normal 
lead. As expected G NN (9) is a monotonic function of 
9, whereas G NS (9) possesses an extremum at 9 C ~ tt/8. 
To understand why the extremum is a maximum, recall 
that for 9 = or 9 = tt, when spin is conserved, current 
flows when a right-going (spin a) electron passes through 
the multilayer, Andreev reflects as a left-going (spin —a) 
hole, which retraverses the multilayer. A M-layer whose 
moment is aligned with the spin of the incident electron 
is anti-aligned with the spin of the outgoing hole and 
consequently the number of aligned and anti-aligned M- 
layers encountered by a given quasi-particle is the same 
for both 9 — and 9 = tt (only the order differs). When 
the elastic mean free path is comparable with the to- 
tal multilayer length, the resistance of traversed layers 
add in series and therefore, apart from small differences 
due to interference effects Jl|], G NS (0) ~ G NS (ir). Fur- 
thermore, since a quasi-particle must necessarily traverse 
regions in which it is a minority spin, both G (0) and 
G NS (ir) are low-conductance states. In contrast, as 9 
increases from zero, this conductance bottleneck is re- 
moved, because an Andreev reflected minority hole can 
spin-convert to a majority hole, thereby avoiding anti- 
aligned moments on its return journey. Of course this 
initial increase in G NS (9) is eventually overcome by the 
usual GMR effect which decreases G NS (9) as 9 — > 7r, 
thereby producing an overall maximum. 
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In the absence of disorder, the nature of the extremum 
is determined by interface scattering and band struc- 
ture. To illustrate this consider a clean multilayer which 
is perfectly periodic and therefore the variation of the 
conductance with 9 arises from tuning of the ballistic 
spin-filtering by the structure. Fig. |^ shows the conduc- 
tance divided by the number of open channels in the left- 
hand-side normal lead. As expected, G NN (6) is a mono- 
tonic function of 9, whereas G NS (6) exhibits a minimum 
around n/2 and then increases. (In this case, transla- 
tional invariance in the transverse direction allowed us 
to use a full ab initio, spd Hamiltonian to obtain the re- 
sults of Fig. ||.) In the NN case, the dependence of the 
multilayer resistance on 9 predicted by our model is in 
good agreement with experiment jL(| . In Ref. the ra- 
tio between the resistance at a given 9 and the resistance 
with AP alignment has been found to fit the following 
function: 

= l-acos 2 (0/2) + 6cos 4 (0/2) (5) 

R(ir) 

where a and b are fitting constants. In Fig. [| we show 
the plot of such a ratio for the disordered multilayer con- 
sidered above, along with the best fit to function (||). In 
addition we also checked that this ratio cannot be fitted 
with the same accuracy assuming a pure dependence on 
cos 2 (6»/2) {i.e. with b = 0). For G NS (9) however, no 
such analytic results currently exist. 

Let us now turn attention to the effect of SO coupling. 
Figs. |^a and |5|b show the conductances as a function of 
the SO strength V\ for, respectively, the NN and the NS 
case. In the NN case (Fig. |a) Gp N decreases as Vi 
increases and eventually joins the curve for G^p- This 
can be understood in terms of the heuristic model pre- 
sented in Ref. 0, because, as the SO strength increases, 
the average length required for a spin to flip (spin relax- 
ation length A s /) gets shorter. Therefore in the P align- 
ment, an injected majority electron travels through the 
multilayer for a length A s / before being scattered into a 
minority spin, thereby producing a decrease in the con- 
ductance. This suggests that the value of V\ for which 
Gp N ~ G^p corresponds to a spin relaxation length A s / 
close to the period lp of the multilayer. We have carried 
out a range of simulations which show that this value of 
V\ does not depend on the overall length of the multi- 
layer, but decreases with increasing lp. As expected, the 
conductance with AP alignment does not change signifi- 
cantly with Vi . 

In the NS case (Fig. ||b) the conductance in the P 
aligned state rapidly increases with V\, reaching a max- 
imum and thereafter decreases, eventually joining the 
curve for the AP configuration. Clearly the enhancement 
in Gp S is produced by the onset of spin-flip scattering. 
The abrupt increase is understandable, since even a rel- 
atively small probability for spin flipping opens a highly 
conductive "channel" if the spin-flip events take place in 



the vicinity of the interface. As one can see in the insert 
of Fig. Ill for larger values of V\ , the conductance Gp S 
joins Gp and together they decrease thereafter. As in 
the normal case G^p depends weakly on V\. The value 
Vi ~ 0.08, at which G NS is maximum, corresponds to a 
spin relaxation length close to the total length of the mul- 
tilayer and, as expected, separate simulations show that 
this value of V\ decreases with increasing total length. 
Similarly the value of V\ at which the GMR ratio (of 
Fig. [I]) vanishes corresponds to a spin-relaxation length 
of the order the bilayer thickness lp and is independent 
of the total length of the system. 

In conclusion, we have demonstrated that spin-mixing 
plays a crucial role in determining both the qualitative 
and quantitative features of GMR in magnetic multilay- 
ers with a S-contact. In contrast with the normal case, 
where spin-mixing suppresses GMR, we find that the 
GMR ratio can be dramatically enhanced by the presence 
of spin-orbit interactions and/or non-collinear magnetic 
moments. In experiments carried out to-date, the pres- 
ence of large spin-orbit scattering presumably masks 
the mechanism shown in Fig. |[ which is predicted to 
be a generic feature in the absence of other spin-mixing 
processes. This suggests that lighter metals and super- 
conductors would be more appropriate for observing the 
new extrema predicted in this Letter. Finally we note 
that for the future it would be of interest to examine 
spin-mixing in non-diffusive NS structures such as clean 
spin- valves jt!|, where the GMR ratio can be non-zero 
or negative, even in the absence of spin- flip processes. 




.1001 ■ 1 < 1 ■ 1 < 1 

0.00 0.05 0.10 0.15 0.20 

V 1 (eV) (SO int. strength) 
FIG. 1. Conventional GMR ratio as a function of the 
SO interaction strength for NN and NS cases. Results corre- 
spond to disordered multilayers (W = 0.6 eV) comprising 40 
bilayers with Im = 15 and In — 8. Samples are formed by 
repeating (3x3) disordered unit cells in the transverse plane 
and summing over 25 k-points in the 2-dimensional Brillouin 
zone. The points are an average over 20 disorder realizations 
and the error bars represent the standard deviations from the 
mean. In the insert, comparison between the conductances in 
the P alignment for the NN and NS cases as functions of the 
SO interaction strength. 
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FIG. 2. ^-dependent GMR ratio for NN and NS cases 
in the absence of SO coupling. Results correspond to disor- 
dered multilayers (W = 0.6 eV) of 22 bilayers with Im = 30, 
In = 16, considering a (3x3) unit cell in the transverse plane, 
and a sum over 25 k-points in the 2-dimensional Brillouin 
zone. The points are the average over 50 realizations of disor- 
der. In the insert, ^-dependent conductance for NN and NS 
cases in the absence of SO coupling. 




FIG. 3. ^-dependent conductance for NN and NS cases in 
the absence of SO coupling for a clean multilayer. The multi- 
layer is modelled by a material-specific spd-band Hamiltonian 
(see Ref. [1]), with Co as M-material, Cu as N-material and 
Pb as S-material. Im = 7 and In = 10, considering a (lxl) 
unit cell in the transverse plane, summing over about 5000 
k-points in the 2-dimensional Brillouin zone. 
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FIG. 4. Plot of the ratio along with the best fit 

P 

to the function (0). The value of the fitting parameters are: 
a — 0.407, b — —0.238. Results correspond to a disordered 
multilayer with the same parameters as in Fig. ^| 
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FIG. 5. Conductances for P and AP alignments as func- 
tions of the SO interaction strength for the NN case (a) and 
for the NS case (b). Results correspond to a disordered mul- 
tilayer with the same parameters as in Fig. [lj 
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